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Abstract 

We study a spherically symmetric solution to the Einstein equa¬ 
tions in which the source, which we call an OR-star (orbiting radiation 
star), is a compact object consisting of freely falling null particles. The 
solution avoids quantum scale regimes and hence neither relies upon 
nor ignores the interaction of quantum mechanics and gravitation. The 
OR-star spacetime exhibits a deep gravitational well yet remains sin¬ 
gularity free. In fact, it is geometrically flat in the vicinity of the 
origin, with the flat region being of any desirable scale. The solution 
is observationally distinct from a black hole because a photon from 
infinity aimed at an OR-star escapes to infinity with a time delay. 


1 The problem 

Can photons form a star? Or, more generally, can a collection of null particles 
that interact only via their mutual gravitation form a bound object? 

Consider first a simpler setting: a swarm of Newtonian particles occupy¬ 
ing spherically and temporally symmetric orbits. The classical gravity felt 
by a particle at radius r is due only to whatever mass is inside radius r. 
Therefore: 
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1. If the swarm has an inner edge r inner , then a particle at r inner must be 
moving tangent to that edge since the particle experiences no gravity 
and continues to move (instantaneously) along its straight tangent line 
to a larger radius. 

2. As the particle’s radial position increases, the amount of mass enclosed 
at smaller r values also increases, causing the particle’s path to curve. 
If the mass inside the particle’s radius becomes sufficient, the trajec¬ 
tory eventually curves enough to deflect the particle inwards at some 
maximum radius r oute r, repeating this cycle. 

This intuition carries over to a relativistic version with a swarm of null 
particles, each of which follows a bound null geodesic of the spacetime gener¬ 
ated by the swarm itself. We call the resulting object an “Orbiting Radiation 


Geodesic of free and trapped photons 



Figure 1: A plot showing photons (null test particles) in the spacetime of a 
high-flux OR-star with a center at the origin. The dotted circle represents 
both the (almost degenerate) inner and outer radius of a high flux OR.-star. 
Photons start at the edge of the flat inner radius with various orientations. 
The blue lines show scattered photons that escape to infinity. The green 
lines show trapped photon trajectories that repeatedly cross the the OR-star 
surface in both directions. 
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star” (or OR-star). In this paper, we demonstrate the existence of OR-star 
solutions to the Einstein equations and discuss their properties. Figure [l] 
probes an OR-star solution with test-particle photons to illustrate the geom¬ 
etry of the solution. 

OR-stars differ quantitatively from but share several qualitative features 
with the Newtonian scenario above. The relativistic null swarm has both 
an inner radius ri nner and an outer radius r outer . The inner radius can be of 
arbitrary real dimension, with a flat spacetime geometry when r < r - mn pr - At 
the outer radius, which can be significantly larger than the inner radius, the 
solution matches onto an exterior Schwarzschild geometry. In the OR-star 
region r inner < r < r outer , the null particles follow planar eccentric bound 
null geodesics arranged in a spherically symmetric configuration. As with 
most eccentric bound relativistic orbits, these orbits are not simple ellipses, 
although they do have well-defined pericenters and apocenters at ri nner and 
r outer, respectively. 

To preview some results, the minimum outer radius of an OR-star is 9/8 
of its Schwarzschild radius. And, unlike with Schwarzschild black holes, the 
interior geometry of an OR-star can be probed by observing the deflection 
pattern and propagation delay of test particle photons sent toward the OR- 
star. 

1.1 Contrast with other solutions 

To elucidate some features of OR-stars, it is helpful to contrast them with 
two other families of solutions to the Einstein equations that exhibit compact 
objects with deep gravitational wells. 

1.1.1 Black holes and similar 

The conventional Schwarzschild black hole solution to the Einstein equations 
has a singularity at the origin and an event horizon at the Schwarzschild 
radius. Both of these elements ignore the potential for quantum phenomena 
to modify the solution in the neighborhood of these spacetime regions. 

Gr avast ars 0H2] are black hole alternatives that rely on Planck length 
scale effects to resist gravitational collapse. Boson stars EE] rely on low- 
mass bosons to form compact objects. For particularly large objects, the 
wavelength of the boson constituents must be particularly large. Related 
to boson stars are Geons [13]. Spherically symmetric Geons are known to 
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be unstable, since a lower energy configuration exists with two concentric 
counter-rotating rings PI- The timescale of this instability, however, is 
unclear and almost certainly depends inversely on the size of the object. All 
of these solutions can be of arbitrary scale, similar to an OR-star, but either 
rely on quantum effects or ignore them in regimes where quantum effects 
plausibly matter. 

In contrast, OR-stars do not rely on quantum effects or ignore them in 
regimes where they matter. Its constituents can be any null particle that 
interacts on large scales only via gravitational effects. These particles can 
have any energy distribution so long as the energy of any single constituent 
particle is a negligible fraction of the mass-energy of the entire ensemble. The 
OR-star has a well-defined continuum limit as an anisotropic but otherwise 


ideal fluid, as discussed in Section 3.2 


1.1.2 Other null particle solutions 

Many null particle solutions to the Einstein equations HI SI El 0 EU have 
been studied. Amongst all these solutions, the scenario of colliding radially 
ingoing and outgoing null particle streams [1] is a limiting case of the OR-star 
solution with some substantial technical similarities discussed in Section 17.21 
This solution, however, requires a negative mass at the origin. Although 
OR-star solutions do also admit a negative central mass and can use that 
feature to enhance their similarity to Schwarzschild black holes, OR-star 
solutions do not require exotic negative mass anywhere and in fact can join 
onto a flat spacetime interior solution. See section 0 for details. 


1.2 Outline 

The rest of this paper is organized as follows. Section [2] lays out preliminaries 
like coordinate choices and the form of the Einstein tensor for our problem. 
In Section [3j we derive the stress-energy tensor for an OR-star. In Section [IJ 
we reduce the Einstein equations to a set of coupled ordinary differential 
equations for the OR-star geometry and for the geodesics of the null particles 
that comprise it, deferring several computational details to appendices. 

In section [5j we discuss various approaches for numerical integration of 
those equations. We have tried them all and have verified that their re¬ 
sults agree within numerical error (we have made all our integration codes 
available). 
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In Section [6j we discuss the structure of the solution from several view¬ 
points. Finally, in Section [7J we discuss variations on the solution family 
which can be easily handled and their implications. 


2 The Einstein tensor 


The solution we seek is spherically symmetric and static, so we use radial 
Schwarzschild coordinates with the following notation for the metric: 

ds 2 = —A(r)dt 2 + B(r)dr 2 + r 2 (d0 2 + sin 2 (6*)d0 2 ) . (1) 


The metric functions A(r) and B(r ) should match onto their Schwarzschild 
values at the outer edge r = r outer of the OR-star. 

Because an OR-star is made up of null particles, its total stress-energy 
tensor will be traceless (see equations (21) (23)). The Einstein tensor must 
therefore also be traceless, which implies that the Ricci scalar vanishes iden¬ 
tically. As a result, just as in vacuum spacetimes, the Einstein tensor in an 
OR-star spacetime is simply the Ricci tensor, which in these coordinate^] 
has nonvanishing components 


Rtt 

Ree 

R(j>cf) 


A" A' A' f A' B'\ 

2B + VB~ 1b\A + ~B ) 
_A" B' A' (A! B'\ 

~2A + VB + lA\A + ~B ) 
] 1 _ r (A! _ B'\ 

~B~ 2B\A~ ~BJ 
sin 2 {6) Ree 


(2a) 

(2b) 

(2c) 

(2d) 


3 The stress-energy tensor 

The particles comprising an OR-star follow null geodesics of the spacetime 
generated by those same null particles. A description of those null geodesics 
allows us to deduce the form of the stress-energy tensor of the entire OR star. 
When all those geodesics have the same magnitude impact parameter \b\ and 

1 See, for example, m- Page 363 of this textbook derives these equations using the 

same notation we use here. 
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are distributed in a spherically symmetric way, the resulting OR-star can 
be interpreted as an anisotropic fluid with a distinct radial and tangential 
pressure at every point. 


3.1 Null geodesics of OR-star particles 

We derive the equations of motion along the fiducial null geodesic of a single 
OR.-star particle. The trajectory and 4-velocity of any other particle can be 
derived from this fiducial geodesic by some combination of rotation within 
the orbital plane, rotation of the orbital plane, and a translation of orbital 
phase. 

Spherical symmetry lets us choose the equatorial plane 6 = 7t/2 to be the 
orbital plane of the fiducial geodesic without loss of generality. is then a 
Killing vector, so a particle with 4-momentum p M has a conserved angular 
momentum 

£ = (3) 

that we allow to take any sign. Because it is static, an OR-star spacetime also 
has a timelike Killing vector dt that yields a conserved “energy at infinity” 

. (4) 


We define the impact parameter of each particle to be 



E°° 


(5) 


and allow it to take any sign (matching the sign of L and 0). Because two null 
particles with the same b but different energies still follow the same geodesic, 
it is useful to define an affine parameter A along null geodesics 


dx M 
~d\ 


= = 



( 6 ) 


that absorbs this energy factor so that it need not appear explicitly in the 
equations of motion. 

Equation (jdj) then yields the equation for t, 

E°° = -(-A(r))Eri—>i=-E, (7) 
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while (J3I) yields 0, 


L = r 2 E' 


2 rpoo, 


( 8 ) 


For the fiducial geodesic, 0 = 0 by our choice of equatorial plane. Finally, 
the radial equation of motion follows from the fact that the 4-velocity of 
OR-star particles is null. Using ( 7j) and (JI]) to eliminate i and </>, we get 


7i\ 


0 = g lw u ,l u v = — At" + Br + r~9 + r sin (0 = — )</> 

1 R-2 ^ 

= —r + Br 2 + — 


A 


r 2 = 


r = 


AB 


Ab 2 


Wab 


1 - 


Ab 2 \ 


1/2 


■ 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


Thus the 4-velocity of our hducial null geodesic is 




1/A 


i / 1 /i Ab 2 

^ ' 1 AB V r 2 

0 

b 


1 

A 


| / A. -1 A6 2 

^ A ' B Y 1 r 2 

o 

Ab 

r 2 sin(0) 


(13) 


sin(0)r 2 

More generally, if the spatial velocity of an OR.-star particle at r makes an 
angle 7 when projected onto the sphere of radius r, then its 4-velocity is 
(where 7 = 0 would be latitudinal and 7 = n/2 would be longitudinal) 



1 


\ 1 1 

9 , 0; 7) = j 

sta(7)f 

COs(7) ■ jw; 2 

L ^ 1 / sm( 0 )r^ J 

1 

“ A 

V 


(14) 


with v denoting the ordinary spatial velocity 3-vector in the standard Schwarzschild 
coordinate basis. 


3.2 OR-stars as anisotropic fluids 

The stress-energy tensor for the entire OR-star should be the sum of the 
individual stress-energy tensors of all the constituent null particles. T^ u for 
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a single null particle with 4-velocity u is essentially that of a null dust 


T ~ pu <E> u 


(15) 


with delta functions as appropriate in the energy density p. 

Consider a local orthonormal frame held with basis vectors and basis 
one-forms 


ex = 


Cr — 


d t 

d r 


ea = 


e x = 


\/W) 

1 

7W) 


-d 0 

r 

r sm 6 


de l = —\J A(r)dt 
de r = y/B(r)dr 
de° = rdO 

de^ — r sin 9d(f). (16) 


Since all the particles have the same impact parameter magnitude \b\, the 
distribution in a local orthonormal frame is anisotropic. The projection v 
of v into the spatial plane normal to ef is isotropic within that plane, so that 


<(«T) = <(«*) 2 > = • 


But, as we can see from (14) and (16), the angle between v r and ir 
with r, so that in general, 


(17) 

varies 


<(«') 2 > = 

Thus, the radial pressure 

and the tangential pressure 

Pe = 


i - ((^) 2 > + \(^T) ■ 

Pr = 


P* = 7,P(( V± ) 2 ) = P ± 


(18) 

(19) 


( 20 ) 


are unequal in general. 

Taking the continuum limit, an OR-star is an anistropic fluid whose stress- 
energy tensor in the local orthonormal rest frame of a fluid element is 


T^° — diag(p, P r , P±, P ±), 


( 21 ) 
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or, in manifestly covariant form (see, e.g., |9] or 0 ), 


=(p + P ± )UW + P^ V + [P r ~ 


(22) 


Even though the fluid consists of null particles, each resulting fluid ele¬ 
ment has a timelike 4-velocity tA and a well-defined rest frame. However, 
because the OR-star consists of null particles, we still expect this to be 
traceless. And it is - the speed of light is unity in all local Lorentz frames, 


so from (19) and (20, we get 


p r + 2P ± = p(((vy) + ((iPr) = p. 


(23) 


3.3 Explicit expression for T jJV 

To conclude this section, we derive an explicit expression for the stress-energy 
tensor (J22| in terms of r, the (still unknown) metric coefficient functions A(r) 
and B(r). and a small number of parameters that characterize an OR-star. 

Let n + (r) and n_(r) denote the number densities at r (as measured in 
frame (16)) of null particles with v r > 0 (outward bound) and v r < 0 (inward 
bound), respectively. n + (r) must equal n_(r) at every r - if, say, n + were 
larger, then when the radial velocities of all the particles changed sign after 
half a radial period of the null geodesics, n_ would be larger, counter to the 
assumption of staticity. 


As measured in frame (16), the radial component of the spatial velocity v 


for null particles moving radially outward ( v r > 0 )across the spherical surface 
at r is (using (Jl4|) and (16)) 


v r = df/dt = 


\[B dr , . .. 

\fA dt 

(24) 

Vb y/A 1 Ab 2 

(25) 

y/A \[B V r 2 

1 Ab 2 

\ r 2 ' 

(26) 


In the local orthonormal frame (16), these radially outward-bound particles 
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cross the spherical surface at r at a rate 


dN 2 f 

—— = n_|_(r)47rr v 
dt 


= n + (r)47rr 2 \ 1 


Ab 2 

r-2 ’ 


(27) 

(28) 


where n + (r ) denotes the (also locally measured) number density of outward 
flowing particles. Because the system must be time-reversal invariant, we 
must have 


n + {r ) = n_(r) 



(29) 


where n{r ) is the total local particle number density (outbound and inbound). 

In terms of the Schwarzschild time coordinate t, the rate corresponding 
to (27) is time-dilated to 


dN + 

dt 


Va^!± 

dt 

y/An + (r)47rr 2 u r 


'/A-n(r)A'Kr 2 \j 1 — 


Ab 2 


(30) 

(31) 

(32) 


By staticity, this rate must be independent of t. It must also be independent 
of r - after one coordinate radial period T r of the geodesics, every particle 
returns with the same velocity to the spherical surface from which it departed, 
and the “geodesic conveyor belt” has carried every particle in the OR-star 
outward across the spherical surface at r exactly once. Therefore, if the total 
number of particles in the OR-star is N, then 


2n \J A(r)n{r)r~ 


1 - 


A(r)b 2 


N 

T~ ’ 
± r 


and the local particle number density is 

N/T r 


n(r) = 


2vr 


\/.4(r)r 2 0 - 


(33) 


(34) 


Now multiply both sides of (34) by the average particle energy in the local 
frame ( E loc ) = (E°°) /y 7 A{r). On the LHS, we get the local energy density 


p(r ) = n(r) (E loc ) (r) 


(35) 
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that appears in (22). On the RHS, we get 
(E°°)N/T r 1 


2ir vw) - w 

Mqr ,/T r 1 


2vr 


A(r)r 2 


1 - 


A{r)b 2 


(36) 

(37) 


where Mor is the total mass of the OR-star measured by distant observers, 
and Moft/T r is the “energy circulation rate” of particles within the OR-star. 
We absorb all the prefactors into a single constant 

Mok/Tt 


m = 

and define a convenience function 

m(r) = 


2vr 


m 


r 2 y/l — b 2 A(r)/r 2 
The local energy density in the OR-star is then 

in m(r) 


P(r) = 


A(r)r 2 \J 1 


A(r)b 2 A(r) 


(38) 

(39) 

(40) 


Using (40), (19), (20) and (14), the pressures P r (r) and Pj_(r) can also be 


expressed in terms of r, A , B,m(r ) and b. 

Combining those results with the metric (Jl| and equation (22), we ar¬ 


rive at the following expression for the nonvanishing components (in the 
Schwarzschild coordinate basis) of the OR-star stress-energy tensor 


rjntt 

= m(r)/A 2 

(41a) 

rjnrr 

- ^(i-^V) 

(41b) 

rj-\Q6 

= m(r)6 2 /(2r 4 ) 

(41c) 

r-ptjxf) 

= T ee /sm 2 {6). 

(41d) 

3 tensor = g a ^g^T al3 is 


T t t 

= m{r ) 

(42a) 

T rr 

= Ab 2 /r 2 ) 

A 11 

(42b) 

Tee 

= m(r)b 2 / 2 

(42c) 

Tw 

= Tee sin 2 (6 ) ). 

(42d) 
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Note that an OR-star is specified by two parameters in terms of which all 
other parameters are determined: m (the energy circulation rate in the OR- 
star, which we will sometimes refer to simply as “the flux”), and b 2 (the square 
of the impact parameter). The total mass of the OR-star is also a parameter, 
but adjusting it merely rescales any given solution (see section 4.2). 


4 The Einstein equations for an OR-star 

4.1 Derivation of ODEs 

As discussed in Section [2j the Einstein equations for an OR-star become 


R^u = 8 nGT^ , 


(43) 


with the left- and right-hand sides, respectively, defined by (|2]) and (42). As 
we demonstrate in the Appendix, the fact that 


tr T = T>\, = 0 

A 6 


(44) 


helps us reduce (43) to a pair of coupled, nonlinear lst-order ODEs for the 


metric coefficients A(r ) and B(r ) 


rA'/A 

rB'/B 


fhBy/l — b 2 A/ 

A 

fhB 

Ay/1- b 2 A/r^ 


+ B- 1 


— B + 1, 


(45a) 

(45b) 


where a prime denotes differentiation with respect to r. 

To prove the above, we develop an equation for rA'/A and rB'/B which 
holds for all null particle problems. In other words, only under the assump¬ 
tion that T = g aa T aa = 0, we can get some simplifications of the Einstein 
equations. In particular, we eliminate the A" from our equations for R a and 


Lemma 1. If g uv T uv = 0, we can write our ODE in terms ofT’s: 

A'/A = -l6nGBr 3 T ee + 8nGrBAT u + B/r -l/r (46) 

B'/B = 8nGrBAT a - B/r + l/r (47) 
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The proof is in the appendix. 

Proof of equations (45a) and (45b): From equations 41a|41d , and 
(39) we get: 


T tt = m/A 2 
T ee = mb 2 / (2r 4 ) 


Plugging these in to lemma [l] we get: 

A!/A = — IQiiGBr^mb 2 / (2r A ) + 8nGrBAm/A 2 + B/r — 1/r 
B'/ B = 8nGrBAm/A 2 — B/r + 1/r 


which simplifies as: 

A!/A = —87rGBmb 2 /r + 8nGrBm/A + B/r — l/r 
B'/B = 8-KGrBm/A — B/r + 1/r 


Now plugging in m from equation [39] and setting rh = 8irGm we get: 


which is 


A!/A 
B'/B 


„ rB/A — Bb 2 /r 

771 - - 

r 2 ^/1 — b 2 A/r 2 


+ B/r — 1/r 


rhB 

rA\J\ — b 2 A/r 2 


B/r + 1/r 


A'/A 
B'/B 


and simplifies to 


1 — b 2 A/r 2 

m(Br/A)— — j====== + B/r — 1/r 


r 2 sjl — b 2 A/r 2 


rhB 


rA^Jl — b 2 A/r 2 


B/r + 1/r 


A'/A 
B'/B 


rhB^Jl — b 2 A/r 2 

Ar 

rhB 

rAyJl — b 2 A/r 2 


+ B/r — 1/r 
- B/r + 1/r. 


Multiplying by r completes the proof. 


□ 
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4.2 The solution family 

All OR-star solutions have an inner radius and an outer radius. The inner 
radius ri nner is the smallest radius that constituent null particles achieve while 
the outer radius r outer is the largest radius that constituent null particles 
achieve. 

Other key parameters are the flux rh, the impact parameter b, and the 
boundary conditions for A and B. Using symmetries and dependencies, we 
reduce this to a one-dimensional family of solutions. To see this, we first 
describe two invariances of the ODEs. 


1. If A, and B are solutions for (45a) and (45b) with parameters b and 
rh, then for any constant k, kA and B are solutions for b/\fk and krh. 


A{r),B{r ) satisfy (|45a|)-(|45b|) ==> kA(r), B(r) satisfy (|45a|)-(|45b|) 

(4S) 

We use this property in our simulations to set A(ri nner ) = 1. After the 
fact, this symmetry can be used to set A(oo) = 1 as per convention. 

2. If A(r) and B(r) are solutions ( |45a ) and (45b) with parameters b and 
rh, then for any constant k, A{kr) and B{kr ) are solutions for kb and 
krh. 


A(r),B(r ) satisfy (45a)-(45b) =>- A{kr), B{kr) satisfy (45a)-(45b) 

(49) 

This scale invariance implies we can set r \ nn Pr = 1 in our simulations. 


A few other simplifications are immediate. 

1. We have B{r- mneT ) = 1 since we are seeking a solution which is flat 
inside the inner radius. As long as the mass distribution is compact, 
this implies B{oo) = 1 as well. 

2. Since A(r inner ) = 1 and r inner = 1, we must have 6 = 1. 

3. Given A fy nn pr ) = 1, R(r inner ) = 1, r inner = 1, and the flux rh, the value 
of v outer is determined. 


The remaining free parameter is the flux m which indexes the family of 
solutions. 
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Semi-analytic solutions 

There are three special cases which can be approximated analytically. The 
first is low flux. Since low flux implies low curvature, photons starting at the 
inner radius move almost directly away from the center at large r. So for 
large r, it is basically a counter flux of inward photons and outward photons. 
Hence the solution of Gergely [3] is a good approximation. We can define fho 
as the maximum flux for which a closed OR-star is not formed. 

A second analytic solution arises in the close-to-critical case. Here there 
is just barely enough intensity to close the “orbits.” In this case, the photons 
end up spiraling to and away from the outer radius with their radius changing 
infinitesimally, limiting to a perfect circle. The following theorem applies 
with proof in the Appendix. 

Theorem 1. If b = 1 and m close to m 0 , then B(r ) ~ 3 and A(r) ~ r 2 * is a 
solution to our equations. In particular, 



1 — 9 rh 2 r 4 /4 + 0(m 6 r 12 
3 + 9m 2 r _4 /2 + 0(m 4 r~ 8 ) 


(50) 

(51) 


B 


A third analytic solution arises in the high-flux case. If the flux is high 
enough, then the inner and outer radius are almost the same. Hence there 
is basically no dependence on r in the ODEs and they can be solved exactly. 
The following theorem applies with proof in the Appendix. 

Theorem 2. For high flux (i.e. m ~ oo), the maximum achievable value for 
B is B & 9, with r inner ~ r outer ~ 2.25. 

From these two theorems, we can catalog OR.-stars by the maximal value 
that B ever obtains. If B max < 3, then the OR-star must be open and the 
photons never return. If 3 < B rnax < 9 then the OR-star is closed with a 
flat interior. There are no OR.-stars with a flat interior spacetimcj^] that have 


B >9 


max 


5 Methods of Integration and Simulation 


Given the partial differential equations from the previous section, there are 
several ways to do an integration (or simulation) to discover the geometry of 


2 We see later in Section 7.2 that B max > 9 is possible with a negative mass singularity 


max 


at the origin. 
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the object. We tried two in order to convince ourselves that the results are 
sound. Both approaches can be easily implemented on a desktop computer 
and run fairly quickly. 

5.1 Integrating Differential Equations 

At r inner and r outer , where the spatial velocity of the null particles is instan¬ 
taneously parallel to the inner and outer spherical surfaces of the OR-star, 
the density diverges to infinity like p{r) ~ 1/^/r — T inn er/nutpr. But this den¬ 
sity still integrates to a finite value, so that the total mass can be well- 
approximated numerically by integrating very close to r inner and r oute r- We 
tested this approach using the deSolve package in R which uses Runge-Kutta 
techniques for integration. 

5.2 Integration from the inner radius to the outer radius 

To confirm that we are not losing essential physics near the inner and outer 
radius, we also tried an alternative numerical approach, using as the inte¬ 
gration variable the affine parameter of an element of null particles traveling 
from r inner to r outer . 

In essence, this approach does integration-by-simulation. Since spherical 
symmetry exists, we can pick any element at the inner radius and follow its 
trajectory to the outer radius, referring to the element’s radius to index a 
differential shell of elements of two sorts: 

1. Elements starting siuniltaneoushj^] on the inner radius with positions 
uniform on the inner sphere and trajectories drawn from uniform tan¬ 
gents to the inner sphere. 

2. Elements of the same mass/energy at the same radius returning from 
a hypothesized outer radius. Conservation principles imply that such 
elements must have the same mass/energy as elements of the first sort 
differing only in the sign of their radial velocity. 

This simulation is self-consistent if the curvature generated creates an 
outer radius where the geodesic has no radial component. At this point, the 
simulation can halt because we have already accounted for the contribution 
of elements returning from the outer radius to the inner radius. 

Simultaneous is defined with respect to an observer at the origin. 
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In more detail, we can start the simulation with r = 1, .4(1) = 1, B(l) = 1 
at the inner radius and a free parameter given by m as discussed in section 4.2| 
Furthermore, the null particle must be tangential to the inner radius, imply¬ 
ing 6=1. 

If we let 6 be the simulation scale, the equations for the forward step then 
become (after some algebra): 


1 . m <— //the mass inside the radius 

2. dm +- //the change in mass 

3. 0 | ln(.A)//the current angle of null particle 

4. d(f) +- S Arn+ ^^ 2 - A //tiie change in angle 

5. dr +- /the change in radius 

This leads to new values according to: 


1. r +- r + dr 

2. m +- m + dm 


4. 0 +- 0 + d(j) 

5. A +- e 2 ^ 


We can record the values of r, A, B and simulate forward to the outer 
radius. With a self-consistent solution for the shell, attaching a Schwarzschild 
geometry onto the outer radius provides a solution everywhere. Then A(r) 
can be rescaled so A(r) = 1 in the asymptotic rest frame. 

We implemented this integration/simulator using C++ doing simple first- 
order integration according to Euler’s method (as above). The results are 
quantitatively stable to high precision and agree with the Runge-Kutta ap¬ 
proach. Our code is available [3]. 
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5.3 Ray Tracing 

In order to test our understanding of the OR-star object, we also created a 
four-dimensional ray-tracer to probe the family of geodesics that test-particle 
photons traverse |^j This process is a simple computation given the geometry 
A(r) and B(r) and conservation of angular momentum. 

To do this efficiently, we recorded the geometry at each discrete radius 
given by the integration/simulation above and used a simple linear interpo¬ 
lation between the two nearest recorded points. 

Searching for the geometry can be done for a point at an arbitrary radius 
by first testing whether the radius is within the shell of the object. If it’s inte¬ 
rior, the geometry is flat. If it’s exterior, then the geometry is Schwarzschild. 
If it’s in the shell, then a binary search suffices, although a hinted location 
reduces this search to amortized constant time in the case of ray tracing. 

6 Simulation results 

We exercised our simulations in various ways to both check correctness and 
understand the implications of OR-stars. 

An important concept is the characteristic geodesic. A characteristic 
geodesic is the path that null particles takes when starting tangential to the 
inner radius at the inner radius. When this is gravitationally closed, the 
characteristic geodesic can be glued together from inner radius to outer to 
inner to outer etc... to get the path of null particles. 

6.1 Varying initial conditions 

Here we exercise initial conditions to get a sense of the solution structure. 

6.1.1 Varying radius 

Figure [2] shows the structure of the characteristic geodesic is invariant to the 
radius, implying that the scale of the structure can be arbitrary, as expected. 

4 These photons are test-particles for observation only, distinct from the null particle 
constituents of the OR-star. Thus, they need a high enough energy to have a small 
wavelength with respect to the scale of the OR-star but a small enough energy not to 
disturb the solution. 
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Characteristic geodesic of null OR-star particles from inner to outer radius 



Figure 2: A plot of the characteristic geodesic from inner to outer radius for 
OR-stars with varying radius but the same flux. This graphically shows (49). 


6.1.2 Low mass/energy flux 

In Figure [3] we see that the structure of the characteristic geodesic as the 
mass/energy flux varies. This result has more structure, because the results 
with 0.25 or 0.5 flux are not closed. 

This result is as expected—low amounts of mass/energy should not be 
able to generate sufficient curvature for geodesics to curl back. 

In Figure [4] we show the time geometry A = g tt and radial geometry 
B — g rr . The difference between our model and colliding null particles [3] is 
that colliding null particles head directly towards or away from the center— 
an impact parameter of zero. Making this works requires a negative mass 
singularity at the origin. We avoid this with null particles having a greater 
impact parameter, generating an inner radius of closest approach. At a few 
multiples of the inner radius from origin, the photons head almost exactly 
away from the origin so the equations are similar to the colliding null particles 
solution. 

For such an open object, the maximal value of B is always less than 3 
namely the B corresponding to the photon sphere. The open object with 
the largest B has null particles orbiting in a photon sphere at each possible 
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Characteristic geodesic of null OR-star particles from inner to outer radius 



x 


Figure 3: The characteristic geodesic for OR-stars with the same inner radius 
but different flux. 

radius. 

6.1.3 The critical mass/energy flux 

In Figure [5] we see variations around the critical point where curvature 
becomes sufficient to create closed geodesics. In this plot, all non-closed 
geodesics reach the graph edge and are valid when the null particles are 
reflected towards the origin. 

Figure [6] shows A, B , and the potential well created by a near-critical 
closed solution. B just barely increases past 3 and the null particles re¬ 
turns towards zero. Near-critical closed solutions are unstable since photons 
significantly off the characteristic curve escape to infinity. 

6.1.4 High mass/energy flux 

The final example in this section is a high flux version (see figure [7]). The 
maximal B = 1/(1 — 2/2.25) = 2.25/.25 = 9 matches the bound in Theorem 
[2] for the maximal mass which doesn’t collapse into a black hole. This is a 
much more stable object as illustrated by the trapped photons of Figure [l] in 
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r 

Figure 4: An “open” OR-star with insufficient null particles for a compact 
object. The upper plot shows A with a yellow background showing the region 
of null particles. The lower plot shows the value of B and the density of the 
mass. 

the introduction. 

6.2 Observational differences from a black hole 

The structure of the object is similar to but observationally different from 
a black hole because geodesics from infinity can probe the geometry at all 
points and return to infinity. In Figure [8] we show the geodesics traced by 
photons from a point source far away. The expected external structure of 
Schwarzschild geometry is observed, but the internal structure is substan¬ 
tially more complicated. It’s easy to observe that the tics of coordinate time 
advance more slowly in the interior from the viewpoint of an observer far 
from the object, as expected. 
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Characteristic geodesic of null OR-star particles from inner to outer radius 



x 


Figure 5: A plot of the characteristic geodesic from inner to outer radius for 
OR-stars with the same inner radius but varying flux near the critical flux. 
Geodesics reaching the edge of the plot do not close. 


In Figure [9] we zoom in to see the internal structure in the flat internal 
region. As expected, the geodesics are all straight lines in this region, al¬ 
though it is interesting to notice that the geodesics all focus on a point. The 
focal point is well within the inner radius and clearly imperfect. 

Considering the global structure again, the object differs from a black 
hole observationally for any geodesic crossing the outer radius. Every photon 
crossing the outer radius is less tangential than the null particles composing 
the OR-star. This monotonicity property is preserved as the photon and null 
particles travel to the inner radius where the photon must therefore enter the 
inner radius since it is not tangential. 

Photons aimed directly at the center have their coordinate time retarded 
by the mass of the object but are otherwise unaltered as observed asymptot¬ 
ically. Hence, the object acts as a “timescope” allowing an observer to record 
what happened on the other side of the object in the past. 

One observational signature of the OR-star for a stationary observer is 
therefore a combination of normal gravitational lensing via Schwarzschild 
geodesics plus a time delayed observation of a point source. When the point 
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Figure 6: This is a near critical OR-star. The upper plot shows A with a 
yellow background showing the region with null particles. The lower plot 
shows the value of B and the density. 
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Figure 7: This is a high flux OR-star. The upper plot shows A with a yellow 
background showing the region with null particles. The lower plot shows the 
value of B and the density. 
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Geodesic of photons from a point source interacting with low flux OR-star 



x 

Geodesic of photons from a point source interacting with high flux OR-star 



Figure 8: Geodesics as traced by photons from a point source far to the 
left interacting with a near-critical OR-star (above) and high flux OR-star 
(below) at the origin. The origin is marked in red, and the structure of the 
OR-star (inner radius, outer radius, and characteristic geodesic) are shown 
with dotted lines. For the high flux OR-star these overlap leaving only a 
circle. The black, red, and light blue tics mark equal separation in time 
as measured by the point source. The red lines are Schwarzschild geodesics 
while the blue lines are not. 


25 



























Geodesic of photons from a point source interacting with low flux OR-star 



Figure 9: A zoom on the same plot as figure [8] showing that the focal point 
is imperfect and the internal space is flat as expected. 


source and the observer are moving relative to the OR-star, it is still possible 
to observe a signature, because there often exists a geodesic between the 
point source and the observer passing through the object. Observation of 
the geodesic passing through the object is difficult because small changes in 
input angle from the source result in large changes in the required angle of 
an observer. In essence, an OR.-star lights up for observers at many angles 
with only a small solid angle of input photons. Since the small angle of input 
photons are shared by many observers, it is relatively dim. 

Hence, a second signature of an OR-star is given by widely dispersed 
observers that record a dim version of the same time-delayed event nearby 
to the object. On such possibility would be given by a nearby star that 
novas. First, observers would see nova via Schwarzschild geodesics, and then 
a dim time-delayed “echo” of the nova event via non-Schwarzschild geodesics 
passing through the OR.-star. 
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7 Solution Variations 


Several variations on the object are straightforward given the above. We 
discuss them here for further consideration. 

7.1 Mass inside the inner radius 

Up to this point we have been using a flat geometry inside the inner shell, 
but other spherically symmetric inner geometries can be handled as long as 
geodesics starting on the inner sphere increase in radius. The only modifi¬ 
cation to the inner solution is an increased amount of time dilation which 
can be determined after the fact as per our existing simulation/solution. The 
calculation of the characteristic geodesic and the geometry between the inner 
and outer shell is only modified by a change to the initial conditions in the 
inside-outside simulation approach. In essence, instead of starting with a 
mass of 0, you start with the mass of the inner solution using the geometry 
that this implies. 

7.2 An OR-star with a (near) event horizon 

If the inner contains flat space, then the smallest outer radius for a mass 1 
object is 2.25 and the maximum value of B is 9. 

But, if we allow more exotic geometries inside the inner radius then the 
outer radius can be close to 2 with B arbitrarily large, similar to the event 
horizon of a black hole. Consider two exotic geometries: an anti-de Sitter 
and a negative mass singularity at the origin. The anti-de Sitter looks like: 

ds 2 = — (l + kr 2 ) dt 2 + (l + kr 2 ) 1 dr 2 + r 2 dkl 2 n _ 2 

and the negative mass of size M 0 looks like the following geometry which has 
a singularity at zero: 

ds 2 = — (1 + Mo/r) dt 2 + (1 + M 0 /r) _1 dr 2 + r 2 dkl 2 n _ 2 

If we use either of these geometries for the interior, we can make the outer 
shell come arbitrarily close to 2 M. Handling this is straightforward—we can 
simply set B to have some smaller value for the initial condition. 

When the outer sphere and the inner sphere are both close to 2 M, we 
have a thin shell of null particles dividing two vacuum regions similar to an 
anti-de Sitter Gravastar [H]. 
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In an OR-star, the null particles eventually turn around and so we have 
a Schwarzschild exterior. We can make this turn around point occur as close 
to 2 M as is desired by placing a negative mass inside of the inner shell. This 
negative mass may be a singularity at the origin, consist of negative mass 
with an internal volume, or fill the volume of the inner sphere similar to an 
Anti-De Sitter space. 

7.3 A large inner radius 

There is no upper bound to the scale of the inner radius, so the only practical 
upper bound is available mass-energy. A natural question therefore is: What 
would an observer inside the inner radius see? 

1. The internal region has a bounded dimension with a definite border 
where the null particles have an integrable singularity. Interaction with 
the null particles at the border may be more or less easy to observe de¬ 
pending on the form of the null particles, but observations from beyond 
the border are certainly distinguishable due to blueshift/redshift effects 
from A(r) varying. 

2. Due to the internal focal point which an OR.-star creates, an internal 
observer could select different internal locations to gain a naturally 
magnified view of locations in the external universe. 

3. Since the internal A(r) is necessarily different from an external A(r), 
all observations of external events are necessarily blueshifted relative 
to internally generated observations. 

7.4 Varying particle energies 

The energy of the constituent null particles was never explicitly used in our 
equations. Consequently, null particles with varying energies could be used 
to construct an OR.-star so long as the individual null particle energies are 
a negligible fraction of the total mass Mor of the OR-star. The negligible 
fraction constraint is imposed by the assumption that null particles have 
uniform directions in the tangent plane at the inner radius. 
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8 Conclusion and Future Work 


The Orbiting Radiation star is a compact star based on null particles capa¬ 
ble of explaining large-mass gravitational wells of any scale. The constituent 
null particles could be something exotic or something more common, like pho¬ 
tons. The OR-star is observationally distinct from most other compact stars 
because test-particle photons entering the OR-star from infinity eventually 
emerge. 

One element not addressed here is formation. Could an OR-star be cre¬ 
ated by a natural or plausible designed process? It is plausible that variations 
of the OR-star can consist of null particles with a variety of impact parame¬ 
ters rather than a single impact parameter since existing solutions can trap 
photons with other impact parameters as in Figure [lj As a consequence, a 
plausible formation process may be somewhat “dirty”, involving null parti¬ 
cles with imprecise impact parameters. However, modeling the dynamics of 
formation and general stability requires substantially more thought. 


References 

[1] J Bicak and P Hajicek, “Canonical theory of spherically symmetric space- 
times with cross-streaming null dusts”, Physical Review D, 2003. 

[2] M Colpi, SL Shapiro, and I Wasserman, “Boson stars: Gravitational equi¬ 
libria of self-interacting scalar fields”, Physical review letters, 1986. 

[3] D Foster and J Langford, OR.-star code, https://github.com/ 
J ohnLangf ord/OR-st ar, 

[4] LA Gergely, “Spherically symmetric static solution for colliding null dust”, 
http://arxiv.org/abs/gr-qc/9809024, 1998. 

[5] H-J Schmidt and F Homann, “Photon Stars”, General Relativity and 
Gravitation 32, page 919, 2000. 

[6] P Jetzer, “Boson stars”, Physics Reports, 1992. 

[7] J Kijowski and E Czuchry, “Dynamics of a self gravitating light-like shell 
with spherical symmetry”, Classical and Quantum Gravity, 2010. 


29 





[8] PO Mazur and E Mottola, “Gravitational condensate stars: an alternative 
to black holes”, http://arxiv.org/pdf/gr-qc/0310107.pdf. 

[9] R Sharma and SD Maharaj, MNR.AS 375, 1265-1268, 2007. 

[10] R Tolman, “Relativity, Thermodynamics, and Cosmology”, Oxford Uni¬ 
versity Press, Clarendon Press, 1934. 

[11] PC Vaidya, “Radiating Spherical Star”, Nature (London) 171, page 260, 


1953. 


[12] M Viser and DL Wiltshire, “Stable Gravastars-an alternative to black 
holes?”, Classical and Quantum Gravity, 2004. 

[13] JA Wheeler, “Geons”, Phys. Rev. 97, 511, January 1955. 

[14] A Zee (2013). Einstein gravity in a nutshell. Princeton University Press. 

Appendix 

proof of Lemma [T] 

First we prove the following: 

Lemma 2. For null dust (i.e. assuming R = g vv R uv = 0) we can define 
S = A'/A + B’/B and D = B’/B — A'/A as the sum and difference of 
log (B) 1 and log(A)', and get the following simplifications: 


S = rR rr + rBR u /A 
D = {2BR ee -2B + 2)/r 
A' = A(S-D)/2 
B' = B(S + D)/2 


(52) 

(53) 

(54) 

(55) 


Proof: Starting with: 



2 BR ee = 2B-2 + rD 
2 ( k BRqq + 1 — B)/r = D 
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To derive our equation S we start with R tt and R rr : 

A" A' A' (A’ B'\ 

tt ~ 2B + 7B~ 4B\A + ~B ) 

A" B' A' /A' B'\ 

rr ~ ~2A + VB + AA\A + ~B ) 

Substituting in S and multiplying by 4 B and 4 A respectively we get: 

4 BRa = 2A" + 4A'/r - A'S 
4 AR rr = —2 A" + 4AB'/(rB) + A'S 

Adding these together we get: 

4 (AR^+ B Rtt) = ±A!/r + 4AB'/(rB) 
r(AR„. + BR tt )/A = A'/A + B'/B 
r(AR rr + BR tt )/A = S 

□ 

Proof of lemma [l| Define S = A'/A + B'/B and D = B'/B — A'/A as in 
lemma [2| and using what we learned there: 

S = rR rr + rBR tt /A 
= r B (R rr / B + Rtt / ’A) 

= rB(R rr B + R tt A) 

= 8nGrB(T rr B + T tt A) 


Using the fact that T is traceless, so —AT tt + BT rr + 2r 2 T ee = 0 we get: 


S = 8TrGrB(2AT tt - 2r 2 T ee ) 


Continuing with D\ 
D 


S + D 
(S + D)/2 
(S — D)/2 


(2 BR ee - 2B + 2)/r 

(2 BR 6e r A -2 B + 2)/r 

(16vr GBr 4 T ee -2 B + 2)/r 

16nGrBAT tt — 2 B/r + 2/r 

8’KGrBAT tt — B/r + l/r 

-16v TGBr 3 T" + 8nGrBAT tt + B/r - l/r 


□ 
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Proof of Theorem CD 

This section proves theorem [l] which derives the limiting solution to our 
equations when 6 = 1 and rh — > m q. Define, 

2 

F = -Jl-A/r 2 (56) 

m 

Lemma 3. Our ODEs are: 


3 r 4 + m 2 F 2 — 2 fh 2 FF'r 
rB' 


B{m 2 F + r 4 — fh 2 F 2 ) 
B 2 

(1 — fh 2 F 2 /r 4 )F 


Proof of Lemma [3} Taking b = 1 we have: 


rA'/A 

rB'/B 


rhBy/ 1 — A/ 

A 

friB 

Ay/1 - A/r 2 


+ B- 


-B + 1 


1 


From our definition of F = 



A 

A' 

rA! 

2 r 2 — 2rh 2 FF'r ~ 1 + 2 m 2 F 2 r^ 2 
2 r 4 — 2 rh 2 FF'r + 2 rh 2 F 2 
2 r 4 — 2 rh 2 FF'r + 2fh 2 F 2 
3 r 4 + fh 2 F 2 — 2fh 2 FF'r 


A/r 2 we get: 
r — m r /r 

2 r — 2 m 2 FF'r~ 2 + 2 rh 2 F 2 r~ 3 
rfiBy / 1 — A/r 2 + (5 — 1)A 
m 2 BF/r 2 + (B — l)(r 2 — fh 2 F 2 /r 2 ) 
m 2 BF + (B- 1 )(r 4 - m 2 F 2 ) 
fh 2 BF + Br A — Bfh 2 F 2 — r 4 + m 2 F 2 
B[fh 2 F + r 4 — fh 2 F 2 ) 


Now working on the second ODE: 


rB'/B 
rB'/B 
rB' 


mB 


Ay/l — A/r 2 

Bfh 


-B + 1 


(r 2 — rh 2 F 2 jr 2 )Fm/r 2 


B + 1 


5 2 


(1 - rh 2 F 2 /r 4 )F 


-B 2 + B 


□ 
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Lemma 4. If lim r 5 = 3, then F and B can be written as polynomials 

—4 

m r . 


Proof of Lemma 0} 

Call the space spanned by polynomials in r~ 4 , V 4 . Define x = r -4 , then 
we are claiming that F — Fq + Fpx + F 2 X 2 + • • • and 5 = 3 + 5i:r + B 2 X 2 + 
■ • •. Technically, what this means is that F and B are analytic in x in the 
neighborhood of x = 0. We won’t take that route. Instead, we show that 
writing 

F = 5 0 + Fir -4 + F 2 r~ 8 H- 


and 

5 = 3 + B\r~ A + 5 2 r -8 + ■ • • 

allows us to solve our ODE’s without introducing terms of the form r _4l+1 , 

4z+2 qj, j,— 4z+3 

First notice that if F e 5 4 then rF' e P 4 . Likewise, if 5 G 5 4 then so is 
r5' G 5 4 . We abuse notation and write any term which is in V 4 as 5 4 . So, 
5 4 + V 4 = V 4 and V 4 * V 4 = Va- So, 


3r 4 + m 2 F 2 — 2 m 2 FF'r 
3 r 4 + fh 2 F 2 — 2 m 2 FF'r 
3r 4 + 5 4 x P 4 — 5 4 x P 4 

V 4 


B(m 2 F + r 4 — m 2 F 2 ) 

Brh 2 F + 5r 4 — Brh 2 F 2 

P 4 x 5 4 + 5r 4 — P 4 x 5 4 x 5 4 

(5 - 3)r 4 


But, since 5 — 3 has a constant term which is 0, we see that (5 — 3)r 4 G V 4 . 
If we now consider our second ODE, 


rB' = 


rB' + B 2 -B = 


5 2 


(1 - fh 2 F 2 /r A )F 
B 2 


B 2 + B 


(1 - fh 2 F 2 /r 4 )F 


5(1 — fh 2 F 2 /r A ){rB' + 5—5) = 5 2 


All terms are in V 4 . 

□ 
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Lemma 5. If F and B are in V 4 and B 0 = 3, then: 


B x = m 2 (F 2 -B 0 (F 0 ^F 2 )) 

B 2 = rh 2 (2FqFi + 8 FqF\ + 2BqFqF\ + B\Fq — BqF\ — B\Ff) 


Fq — B 0 / (. B 0 — 1) 

„ _ 2B 0 B 1 -F 0 (rh 2 F 2 (B 2 -B 0 )- AB 1 +2B 0 B 1 -B 1 ) 

Fl ~ 


Proof of lemma [5} 

F = F 0 + Fir -4 + F 2 r~ 8 H- 

rF' = 0 — 4F 1 r -4 — 8 F 2 r~ 8 + • • • 

B — 3 + Bir~ A + B 2 r~ 8 + • • • 
rB' = 0 — 4F 1 r -4 — 8 B 2 r~ 8 + • • • 

3r 4 + fh 2 F 2 — 2 fh 2 FF'r = B{fh 2 F + r 4 — rh 2 F 2 ) 

Br 4 = 3r 4 + fh 2 F 2 — 2fh 2 FF'r + Brh 2 F 2 — Brh 2 F 
extracting the r 4 terms yields Bo = 3. Extracting the constant terms yields: 

B 1 = m 2 {F 2 - B 0 (F 0 - F 2 )) 

Continuing with the r~ 4 terms: 

B 2 = m 2 {2FoF l + 8F 0 Fi + 2BoFqF\ + B\ F { 2 — BqF\ — B\Fq) 

Proceeding in this fashion we get a series of equations involving higher 
and higher order terms for F*. 
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Now using the other ODE for B we get: 


rB' = 
rFB' = 
rFB' + FB 2 — FB = 


B 2 


(1 - fh 2 F 2 / r 4 )(r F B' + FB 2 - FB) 
F( 1 - rh 2 F 2 / r 4 ) (r B' + B 2 - B) 


F( 1 - m 2 F 2 /r 4 ) 
B 2 

(1 - m 2 F 2 /r 4 ) " 
B 2 

(1 - m 2 F 2 /r 4 ) 

B 2 


-B 2 + B 


FB 2 + F5 


= hr 


Extracting the constant terms out of this we get: 


(F„)(1)(0 + Bl - £t„) = B 2 
B 0 B 0 2 — -^o-Bo = Bq 
FqBq — F 0 = Bq 

Fq = B q / ( B 0 


1) 


Extracting the r 4 terms (This means extracting two linear terms from the 
product along with one r -4 , which can be done in three ways.): 

(FO(1)(0 + S 0 2 -So) 

+(F 0 )(rh 2 F 2 )(0 + B 2 -B 0 ) 

+ (Bo)(l)(—4Bi + 2B 0 Bi — B\) = 2B 0 Bi 

Fi(B 2 — Bq) + Fo(m 2 Fq (Bq — Bq) — 4Bi + 2BqBi — B\) = 2BqB\ 

F\ ( Bq — Bq) = 2BqBi — F 0 (m 2 Fq (Bq — Bq) — 4 Bi + 2B 0 Bi — Bi) 

„ _ 2BqB 1 - F 0 (rh 2 F 2 (B 2 - Bq) - AB X + 2B 0 B 1 - B x ) 

Fl ~ 

□ 


Lemma 6. Solving the equations in lemma\5_\we get: 

Bq = 3 
Fq = 3/2 
Bi = m 2 9/2 

Fi — 0 

B 2 = rh 4 27/8. 
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Proof: We have B 0 = 3, so substituting this in to P 0 : 

F 0 = 3/(3 - 1) = 3/2 
Plugging these both in to D \: 

B 1 = m 2 (( 3/2) 2 - 3((3/2) - (3/2) 2 )) = rh 2 9/2 
Plugging all three into F \: 


6(m 2 9/2) - (3/2)(m 2 (3/2) 2 6 - 4(m 2 9/2) + 6(m 2 9/2) - (m 2 9/2)) 

PI = - 

6 

_ _ 2 27- (3/2)(27/2 - 18 + 27- 9/2) 

1 6 
Pi = 0 


Likewise: 

P 2 = m 2 ((m 2 9/2)(3/2) 2 — (rh 2 9/2)(3/2)) 
P 2 = m 4 27/8 


Proof of Theorem [Q 

We have 

P = 3 + (9/2)m 2 /r 4 + (27/8)?7i 4 /r 8 + 0(r -12 ) 
P = 3/2 + 0(r -8 ) 

A = r 2 — fh 2 F 2 /r 2 
A = r 2 — (9/4)rri 2 /r 2 + 0(r _1 °) 


□ 


□ 


Proof of Theorem [2] 

Before proving this theorem, we first introduce a change of variables: 

z = m 2 {r — 1) (57) 

H = \/2z ~ rh 2 (A — 1) (58) 

We are assuming that the inner turning point is set equal to 1 and that 
.4(1) = 1. The following lemma (whose proof is given later since it is just 
algebra) converts our ODEs to these variables under the assumption that 
m ~ oo. 
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Lemma 7. Define z and H as in (51) and (58). Then after taking the limit 
as fh — > oo we get: 


dH/dz = (3/H - B - B/H)/2 
dB/dz = B 2 /H 


Proof of Lemma [3 

Our original ODEs are (with h = 1): 


rA'/A 

rB'/B 


mBy/l — A/r 2 


A 

fhB 


Ay /1 — A/r 2 


+ B- 1 


B + 1 


Proof of the change of variables: 


dz 

dB/dz 

dH/dz 


m 2 dr 
B' /fh 2 

(1 - A!/2)/H 


This leads to 


A' 

A 


r 

fhry/l — A/r 2 


2 - 2 HdH/dz 
l + (2 z-H 2 )/m 2 
1 + z/rh 2 
rhy/r 2 — A 

m\J (1 + z/rh 2 ) 2 — 1 — (2 z — H 2 )/fh 2 
fhy/ z 2 /fh A + H 2 /in 2 
\J H 2 + z 2 /fh 2 


(59) 

(60) 


Working these into our first ODE: 

rfiBy/l — b 2 A/r 2 


+ BA/r — A/r 


= (B / r 2 ) \/H 2 + z 2 /m 2 + (B — T)A/r 


= (. B/r 2 )y/H 2 + z 2 /m 2 + (B - 1) 


(1 + (2 z-H 2 )/m 2 ) 


A' 

2 - 2HdH/dz 
2 - 2HdH/dz 
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1 + z/m 2 


















Now taking the limit as m —> oo 


2 - 2 HdH/dz 
-2 HdH/dz 
dH/dz 


bVip + b-i 

BH+B-3 
3/(2H)-B/2-B/(2H) 


Now trying the second ODE: 
dB / dz = 

dB/dz = 


Amr\J 1 — b 2 A/r 2 
B 2 

H 2 + z 2 /m 2 

B 2 


dB/dz = - ._ 

(1 + (2z — H 2 )/m 2 )^H 2 + z 2 /m 2 

Now taking the limit as m —> oo 

B 2 

dB/dz = — 

H 


- B 2 /rm 2 + B/rrh 2 
B 2 /rrh 2 + B/rrh 2 

B 2 / ((1 + z/rh 2 )r 


So onr equations are: 


2 HdH/dz 
HdB / dz 


3 - BH - B 
B 2 


Proof of Theorem [2j Taking the ratio of ([59]) to (]60| we get: 


dH 

~dB 


(3 /B 2 - H/B - 1/B)/2 


Which can be solved as for some constant k0 


H = k/y/B-l-3/B 

= -3(l/v / 5-a/3)(l/v A B- 1/a) 

5 a check: 


-kB~ 3/2 + 6/B 2 = 3/B 2 -kB~ 3/2 + 1/B + 3/B 2 -1/B 


2 ) + B/(( 1 + z/m 2 )m 2 ) 


□ 
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where k = a + 3/a. 

In the original variables the turning point condition is b 2 A = r 2 . This 
corresponds to a turning point condition is H — 0 j£] The turning points are 
then: 


r>. 

inner — a 

-^outer = 9/a = 9//ij rmer 

For the flat interior case, we have one root at B = 1 and hence the second 
root is at B = 9 which means k — 4: 

H = A/VB — 1 — 3/B 

Since we know the geometry is Schwarzschild outside r outer , and B = 9 = 
1/(1 — 2/r outer ), we can deduce that r outer = 9/8. By the construction we see 
that r inner ~ v, outer• This completes our claims. □ 

Notice that if we allow negative mass solutions on the inside, any combi¬ 
nation works out for which BiB a = 9. Alternatively, if we consider a black 
hole at the center, then we can have a heavy shell of “orbiting radiation” 
at any distance from 2.25 out to 3.0. This shell provides the extra mass 
necessary to keep them in orbit. 


6 If H = 0, then 2z = rh 2 (A — 1) which is the same as 2(r — 1) = A — 1. Which is 
A = 2r — 1 = r 2 — z 2 /rh 4 . So as m — > oo we see that H = 0 are the turning points. 
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